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Abstract 

Molecular motors convert chemical energy into mechanical work while op- 
erating in an environment dominated by Brownian motion. The aim of this 
paper is to explore the flow of energy between the molecular motors and its 
surroundings, in particular, its efficiency. Based on the Fokker-Planck equation 
with either N or infinite chemical states, we find that the energy efficiency of 
molecular motors, whether the Stokes efficiency or the usual thermodynamic 
efficiency, is strictly less than one, because of the dissipation of the energy in 
both the overdamped surroundings and in the process of the chemical reaction. 
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1 Introduction 

Molecular motors are biogenic force generators acting in the nanometer range 

([12]). Linear motors produce sliding movements along filamentous structures called 
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protein tracks; for example, myosin slides along an actin filament ( [Tf [TT| [27] ) . kinesin 
( [H El [ini EHl EDI) and dynein ([2Sl IS2] ) along a microtubule. Many of molecular 
motors have the ability to rotate; for example, the bacterial fiagellar motor ([HI) and 
-FqFi-ATP synthase ([IS]). This movements and forces generated by the motors play 
essential roles in cellular functions. The linear motor kinesins are widely distributed 
in almost all eukaryotic cells. These motors are adenosine triphosphate (ATP)-driven 
walking machines that move in 8-nm steps toward the plus ends of the microtubule, 
turn over one ATP molecular per step under a range of loads ( [31 El [121 HH] ) • Kinesin 
has two motor domains called heads, each of which includes one ATP- and microtubule 
binding site ( [Il[l5]). This motor steps mainly in the forward direction (to the plus 
end of the microtubule), but occasionally in the backward direction (to its minus end) 

M). 

The study of the energetics of molecular motors is relevant for several reasons. 
Firstly, the efficiency of molecule motors should be high enough to decrease the heat 
dissipation. Otherwise, the temperature of the surrounding would become higher 
and higher since the heat is difficult to be dissipated. Secondly, molecular motors 
are related to fundamental problems of thermodynamics and statistical mechanics. 
For instance, one of the first well-known Brownian motors was the ratchet-and-pawl 
setup studied by Feynman in his Lectures [25j. He calculated the efficiency of such 
an engine since his original purpose was to prove that an automatic demon cannot 
beat the second law of thermodynamics. Thirdly, many models of molecular motors 
proposed in the literature are based on nonequilibrium fluctuations without specifying 
their source. However, the study of the energetics of such models requires a more 
precise formulation, since one has to determine the physical nature of the external 
agent and verify that the motor is consistent with the second law of thermodynamics 

mm)- 

Recently, Hongyun Wang and G. Oster [33j have given the deflnition of the Stokes 
efficiency rjs (see Eq. ( IT7j) ). and shown t]s < 1 through mathematical analysis of the 
Fokker-Planck equation with chemical states. In [23], Hong Qian have introduced 
a nonequilibrium potential function for a motor protein modelled by a rectified Brow- 
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nian motion. Since the conservation of energy can be mathematically established, 
Qian's model provides a comprehensive theory for motor efficiency. In |21j, the en- 
ergetics of forced ratchets, chemical motors, and thermal motors has been discussed. 
The reasonable definitions of the energy efficiency of motors, which are compatible 
with the law of the thermodynamics, have also been reviewed. 

In this research, we will discuss some properties of the energy efficiency of molecu- 
lar motors in more details. Using the description of the Fokker-Planck equation with 

or infinite chemical states, we obtain rjs < l,r]T < 1- It is to say that the energy 
efficiency of molecular motors, either the Stokes efficiency ([31]) or the usual thermo- 
dynamic efficiency ([211 123]), is strictly less than one, and any attempt to design a 
molecular motor that can convert the chemical energy into mechanical energy with 
efficiency one would fail. The movement of molecular motors satisfies the second law 
of the thermodynamics. It is why the Feynman ratchet doesn't work [2H [25] . But 
how is the chemical energy consumed in the movement of the molecular motors? In 
this paper, we'll find that the total energy released by the fuel molecules is consumed 
in three ways: to do mechanical work, to be dissipated into the surroundings, to 
be dissipated in the chemical reaction (see Eq. ( HTi) ). In fact, for motor proteins, 
such as kinesin, myosin and dynein, the heat dissipation in the surroundings is very 
important to keep the body temperature and a well-balanced body temperature is 
essential to the functions of the cell. So, not only mathematical analysis but also 
biophysical demand implies that the energy efficiency of the molecular motor should 
be strictly less than one. For motor proteins, the chemical energy is not only used 
to do mechanical work but also used to keep the body temperature. Of course, high 
efficient molecular motors are also very desirable in mechanical equipment. 

2 Fokker-Planck equation with N chemical states 

In this section, we'll discuss the energy efficiency r]s and t^t using the Fokker- 
Planck equation with chemical states. 

Let Pi{x,t) be the probability density for finding the motor at position x in the 
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i—th chemical state at time t. The evolution of the probability density is governed by 
a set of coupled Fokker-Planck equations that ensure the conservation of probability: 

dpAxA) dJAxA) ^ . ^ . ^ 

^^ + ^^ = /._|(a:,t)-/,^.(x,t) l<^<N (1) 

Here Ji is the probability flux in the spatial direction for the i—th chemical state, 
associated with convective transport and diffusion. Ii+1/2 is the probability flux den- 
sity from state i to state i + 1, is the number of chemical states in one reaction 
cycle. For example, = 4 is usually used in the motility cycle of ATP-driven motor 
proteins ([3ll[3l|). 

In the i—th state, the motor is driven by a periodic potential 0j(a;), and 



Ji{x,t) = pi{x,t)ui{x,t) 



^ dx 



pi{x,t) d(f)i{x) ^dpi{x,t) 



^ dx dx (2) 

Pi{x,t) d(f)i{x) ^ Pi{x,t)fext _^dpi{x,t) 



^ dx ^ dx 



spatial drift flux spatial diffusion flux 

Ii+i{x,t) = pi{x,t)kf — pi+i{x,t)k^_^^ (flux along reaction coordinates) (3) 

where $j = 0j + ksT In pi = (pi — fextx + A;_bT In pi is the enthalpy, ^ is the coefficient 
of viscous drag, D is the free diffusion coefficient, is the Boltzmann constant, T 
is absolute temperature {D,^,kB,T satisfy the Einstein relation ksT = ^D), f^xt is 
the external load {fext < ), and k^' is the transition rate from state i to state i + 1, 
k~^i is the transition rate from state i + 1 to state i. k'f and k~_^^ satisfy 

kt / u±2^i±M\ (4) 



which ensures detailed balance at equilibrium. 

Since the motor operates in a chemical and mechanical cycle, the boundary con- 
ditions for Eq. ([1]) in both the spatial direction and the reaction coordinate are 
periodical: pi{x + L,t) = pi{x,t), pi+N{x,t) = pi{x,t), (pi+N = (pi — A/i, where L is 
the step size of the motor (L = 8.2 nm for motor protein kinesin), A/i > is the 
chemical free energy consumed in one reaction cycle. In physiological conditions, the 
hydrolysis free energy of ATP is A/i ^ 25k bT ([121 EHl 136]). 



At steady state, the dissipation rate Ilj corresponding to slide within the potential 
profile $j (i.e. the dissipation rate in the i—th chemical state) is (P| [2U]) 



n, 



Ji^dx l<i<N 
ox 



(5) 



It can be readily verified that Ilj = i^^Y — ^- Summing ([5]) over i and 

integral by part, we obtain the total spatial dissipation rate 

N ,L N 



i=l i=l 

i=l i=l 



(6) 



In view of Eq. ([T]), at steady state, 

dJi{x) 
dx 



(7) 



Eqs. ([6]) © imply 

i=l "^0 i=l 

i=l ' i=l 



+i{x) - ^,{x))I-^i{x)dx + / Ii{x){^i{x) - ^N+i{x))dx 



i+i{x) - ^i{x))Ii_^i{x)dx 



By Eqs. (|6]) and ([8]), we obtain 

UtL V J^(0)+A/i / Ji ix)dx = V / 

Thanks to Eq. 



p{x) f (9$ 



da; + - $i+i(a;))/i+i (x) 

(9) 



da; 



d ( ^ \ 

— y Ji(x) 1=0 / Ii(x)dx = / Iiix)dx = ■ ■ ■ = / Ij>j,i(x)dx (10) 
J Jo ' Jo ' Jo ^ 



it is to say 



i=l i=l -^0 i=l 



pi{x)ui{x)dx =: y 



(11) 



JO j^]^ ^0 



X =: u 



(12) 



where V is the mean velocity of molecular motors, u is the mean rate of the chemical 
reaction. 

Due to (El) (11 



=Pi{x,t)kf 
=Pi{x,t)kf 
=Pi{x,t)kf 



pi+i{x,t)k 



i+l 



Pi{x,t)k] 



pi{x,t) 



exp 



ksT 



1 — exp 
1 — exp 



so 



- (j)i{x) 

ksT 
^i+i{x) - ^i{x) 
knT 



^i+Ax) -^iix) 



ksT 



^i+iix) - ^i{x) > 0(< 0) =^ exp 

< 0(> 0) =^ (Mx) - > 

Finally, Eqs. © ([II]) ([E]) dHj) mean 

p{x) [d^i^^ 



> l{< 11 



+ {^i{x) - ^i+i{x))^Mx) 



dx 



E 

i=l 



^p{x)ul{x) + {^i{x) - i (x) 



rfa; > 



(13) 



(14) 



(15) 



where / is the dissipation corresponding to slide within the potential profile $j, // is 
the dissipation corresponding to transitions between two chemical states. (ITSl) is the 
second law of the thermodynamics. 

By Eq. ( [I5|) . the thermodynamic efficiency 

fextV 



Vt 



z/A/i 



< 1 



(16) 
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the Stokes efficiency 

i (^J^ pi{x)ui{x)dx 



Vs 



< 



fextV + Z/A/i fextV + Z/A/i 

e/o"Ef=iP.(^)«Ka:)c^x _ Yi:^J^^mYdx (17) 



fextV + Z/A/i /ext^ + Z/A/i 

<1 



In view of ([T] 



{x) - ^i+i{x))I^^i{x)dx = 



:i8) 



^^(<l>,(x) - <l>,+i(x))Ji+i(x) = Vl<i<A^ 0<x<L 

^^<t>i{x) = <^i+i{x) \/l<i<N 0<x<L 

= <l>j(x) Vl<i, j<A^ 0<x<L 

^Afi = 

Therefore, if A/i 7^ 0, then ^fLi J^i^iix) — $i+i(a;))/j_,_i 7^ 0, which imphes 

(see the inequahty (HM and the definitions (fTB]) (IT7|) ) 

r/T < 1 < 1 (19) 

If A/i = 0, then the mean velocity V < (see Eq. (1151) ). It is to say that no useful 
mechanical work can be done by molecular motors if there is no input of the energy. 
If the mean velocity > 0, i.e., molecular motors do useful work against the external 
load fext < 0, then Ajj, 7^ 0, which implies Xlili loi^ii^) ~ ^i+i{x))Ii+i{x)dx 7^ 0, 
consequently r/^ < 1,775 < 1 (see Eqs. fllStllSp ). Therefore, there always exists energy 
dissipation as long as the input energy A/i 7^ 0. 
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3 Fokker-Planck equation with continuous chemi- 
cal coordinates 



In this section, we'll discuss the energy efficiency rjs and riT using the Fokker- 
Planck equation with infinite chemical states. 

Following Refs. fl3[ [22] . we denote y as the coordinate of the internal conforma- 
tional space of the motor protein, and x as its center of mass with respect to the 
linear track. Let p{x, y, t) be the probability density for finding the motor at position 
X and chemical coordinate y at time t, y) be the potential. Due to f|T3l) . in the 
large N limit, 



d 

k-i{x) - I,^ix) = - — 



1 d^x^ 
-^P{x,y,t)^^ 



d_ 

dy 



Dy cM?ix^ 
j^p{x,y,t) 



dy 



(20) 



where D 



lim 



and C = ksT/Dy. Let u{x,y) = v{.x,y) 



•i^^^, in the large limit, Eq. ([T]) arrives 
dp{x,y,t) dJ{x,y,t) dl{x,y,t) 



+ 



+ 



dy 



0<x<L, 0<y <1 



dt dx 
where J = pu, I = pv. 

At the steady state, the mean velocity of the movement is 



(21) 



V =< u >-- 



1 rL 



Jo 



p{x, y)u{x, y)dxdy 



1 rL 



J{x, y)dxdy 



(22) 



the mean rate of the chemical reaction is 

-1 i-L 



V =< V >- 



Jo 



p{x, y)v{x, y)dxdy 



Jo 



1 rL 



Jo 



I{x, y)dxdy 



(23) 



Define 



n 



L rl 



Jo 



d^{x,y) .d^{x,y) 
-J{x,y) — ^ I{x,y)' 



dx 



dy 



I' 



dydx 



(24) 



where /' is the dissipation rate along the spatial coordinate x, II' is the dissipation 
rate along the reaction coordination y. So LI is the total dissipation rate of molecular 
motors. 



It is easy to find 

•1 rL 



at the same time 

-1 rL 



SO 



Since 



~"/o/o J{x,y)^^{x,y)dxdy 
-Jo J{x,y)^{x,y) HZq - Jx{x,y)^{x,y)dx dy 
fextL j'l J (L, y)dy + /^^ J^(x, 
fextL j'l J (L, ?/)c/y - j'l Iy{x, y)^{x, y)dxdy 
fextL Jq J{L,y)dy - I{x,y)^{x,y) {^Zl - I{x,y)^y{x,y)dy 
fextL J (L, y)dy + A/x I{x, l)dx + I{x, y)^y{x, y)dxdy 

n = fextL [ J{L, y)dy + A/i / I{x, l)dx 
Jo Jo 



dx 



d_ 

dx 

d_ 
dy 



J{x,y)dy 



I{x, y)dx 



d d 

J g^-^^^^y^^y ^ ~ J = 

J —I{x,y)dx = -J —J{x,y)dx = 



one knows 



J{x,y)dy = / J{L,y)dy 



therefore 



L rl 



L I J{L,y)dy 
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J{x, y)dydx = V 



^0 
1 pL 



I{x,l)dx = I I I{x,y)dxdy = v 
Finally, (125]) ([271) ([3l]) ([32]) imply 

U = fextV + uAfi = [ [ p {^u^ + Cv^) dxdy > 
Jo Jo 



JO 



(25) 



(26) 



(27) 

(28) 
(29) 



I{x,y)dx= I I{x,l)dx (30) 

(31) 
(32) 

(33) 



this is the second law of the thermodynamics. By (!33|) . the thermodynamic efficiency 



Vt 



fextV _ ^ n 
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< 1 



(34) 



the Stokes efficiency 

^y2 i (^J^ /q p{x, y)u{x, y)dydx^ 



Vs 



fextV + uAjJ fextV + Z/A/i 



lo yWi^, y)dydx (35) 
~ fextV + z/A/i 

_ C Jo Jo p{x, y)v'^{x, y)dydx ^ ^ 

fextV + Z/A/i ~ 

Similar as f|T8|) 

C/ / p{x,y)v^{x,y)dydx = 
Jo Jo 

^^v{x,y) = 0<x<L, 0<y<l 

<^^^{x,y) = ^{x, z) < X < L, 0<y,z<l 

^Afi = <l>{x, 0) - <l>{x, 1) = 

So, if A/i 7^ 0, then ( J^ J^ p{x,y)v'^{x,y)dydx 7^ 0, which means the inequalities 
f|T9l) hold true. 

Moreover, applying the Cauchy-Schwarz inequality, and using the constraint J^ J^ p{x, y)dydx 
1, one obtains 

p{x,y)u^{x,y)dydx> ( / p{x,y)u{x,y)dydx) = 



(36) 



Jo \Jo Jo / 

/ p{x,y)v'^{x,y)dydx > ( / p{x,y)v{x,y)dydx) = z/^ 
^0 V^o ^0 / 



(37) 



so 

fextV + , {^L' + C)V 

Vt = 1— < 1 T = 1 (38) 

z/AjU z/AyU A|UL 

= < ^ ^ 1 _ C .39^ 

here the equahties hold true if and only if u{x,y) = V, v{x,y) = z/. Obviously 
r]T <l, r]s < 1. 

Finally, by Eq. (133|) and the definitions of the energy efficient ris,r]T, we can get 
the following interesting relationship 

^V^ = r]sil-VT)i^Ap (40) 
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= rjrAfi + 775(1 - ?7t) A/i + {1 - t]s){1 - r]T)Afi 

(41) 

From Eq. ( HTil . we can know how the energy is consumed in the movement of molec- 
ular motors: where //tA/i = fextL is the energy used to do useful work against the 
external force fext, (1 — Vt)A^J' is the total energy dissipation, ris{l — riT)Afi is the 
energy dissipation due to the overdamping surroundings, (1 — ris){l — riT)Afi is the 
energy dissipation during the chemical reaction. 

Though the 2D Fokker-Planck equation (12T|) is the large N limit of the equations 
([T]), most of the results in section 3 can't be obtained directly by the limitation of 
the corresponding ones in section 2 (for rj = lim„^+oo ?7n and r/n < 1, it is not always 
correct that 77 < 1). Moreover, through the detailed discussion in section 3, the 
relationship between the parameters in equations ([T]) and the parameters in fl?I]) is 
clarified, and how the energy is consumed in the movement of the molecular motor is 
also clarified easily. The equations ([1]) might be more physical, but it might be more 
convenient to get detailed results using equation (l2Til . 



4 Conclusions 

In this paper, we have given a detailed discussion of the energy efficiency of molec- 
ular motors. We find that the energy efficiency of molecular motors, whether the usual 
thermodynamic efficiency or Stokes efficiency, is strictly less than one. This result 
is sharper than in previous work ([23l El]), where rjs < I and rjx < 1 were proved. 
In each step of a molecular motor, the total energy released by the fuel molecules is 
consumed in three ways: to do mechanical work, to be dissipated in the surround- 
ings, or to be dissipated in the chemical reaction. According to the second law of the 
thermodynamics (!33|) . the faster the molecular motors move, the more the energy is 
dissipated. Theoretically, the highest energy efficiency only can be achieved when the 
movement of the molecular is a quasi-static process, which is more like the Carnot 
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heat engine. 

For the multiple pathways case (0, [ISl [IZ] ), the results in this paper are also 
correct. The only difference is that the energy efficiency (or the energy used to do 
mechanical work, the velocity) in different pathways might be different. Regardless 
of the pathways, the efficiency of molecular motors is always strictly less than one. 

Finally, similar theoretical analysis can be used to describe the backward stepping 
case of motor protein kinesin ([6]). If fuel molecules are resynthesized in the backward 
stepping process, then fext < 0,V < 0, A/i < and fextV + uAji > 0. So the 
energy efficiency of fuel molecules resynthesizing rj = —j^^ is strictly less than one. 
Generally, in each backward step of a molecular motor, the total external energy 
ifextV) is consumed in three ways: to resynthesize fuel molecules (— z/A/i), to be 
dissipated in the surroundings {J^ p^u'^dxdy), or to be dissipated in the chemical 
reaction ( p(v'^dxdy): fextV = — z^A/x + p {^v? + C,v'^) dxdy. Where A/i < 
means fuel molecules are resynthesized, and A/i > means fuel molecules are 
consumed. 
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